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Abstract 

The method for obtaining functional equations, recently proposed by one of the authors 
[1], is applied to one- loop box integrals needed in calculations of radiative corrections to 
heavy-quark production and Bhabha scattering. We present relationships between these 
integrals with different arguments and box integrals with all propagators being massless. 
It turns out that functional equations are rather useful for finding imaginary parts and 
performing analytic continuations of Feynman integrals. For the box master integral needed 
in Bhabha scattering, a new representation in terms of hypergeometric functions admitting 
one-fold integral representation is derived. The hypergeometric representation of a master 
integral for heavy-quark production follows from the functional equation. 
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1 Introduction 



The production of heavy quarks at hadron cohiders has become a very important field of research. 
The very large production rates for both top and bottom quarks at the CERN Large Hadron 
Collider (LHC) will allow for studies of heavy quarks with high precision. The full next-to- 
leading-order (NLO) radiative corrections to the hadroproduction of heavy flavors was completed 
in Ref. [2]. The theoretical NLO predictions suffer from the usual uncertainty resulting from 
the freedom in the choice of renormalization and factorization scales of perturbative QCD. 
To reduce such uncertainties, next-to-next-to-leading order (NNLO) calculations are needed. 
The computation of NNLO corrections is complicated due to great technical difficulties, mainly 
related to the evaluation of two-loop Feynman integrals. To overcome such difficulties, it is of 
great importance to develop new methods and approaches for calculating Feynman integrals. 
The main complications of these calculations are related to the fact that the integrals depend 
on several kinematical variables. As was noted in Ref. [3], the most appropriate methods for 
calculating such integrals may be those based on a different kind of recurrence relations. Such 
methods can be based on the solution of recurrence relations with respect to the exponent of a 
propagator in the integral [4] or on the solution of dimensional recurrences [5, 6]. 

A significant simplification of the computation of Feynman integrals depending on several 
kinematical variables may be achieved by using a new type of relationships between Feynman 
integrals through functional equations with respect to kinematical variables as proposed in 
Ref. [1]. As was shown in Ref. [1], Feynman integrals with several kinematical variables can be 
expressed in terms of integrals with a lesser number of variables, which significantly simplifies 
their evaluation. 

It is the purpose of the present paper to apply the general method for finding functional 
equations [1] to integrals required in calculations of radiative corrections to important physical 
processes and to use those relations for the analytic computation of these integrals. 

Our paper is organized as follows. In Section 2, we give definitions and notations. In Section 
3, we present functional equations for the on-shell master integrals from heavy-quark production 
and Bhabha scattering. In Section 4, a new hypergeometric representation in terms of the Appell 
functions Fi and and the Gauss hypergeometric function for the one-loop box integral 
from Bhabha scattering is presented. In Section 5, we describe how to use functional equations 
to find imaginary parts of the considered integrals. New analytic results for the imaginary parts 
are presented. Using dispersion relations, we write also one-fold integral representations for 
real parts of integrals. In Section 6, we present functional equations for master integrals from 
heavy-quark production with one quark leg off shell. 

2 Definitions and notations 

As was shown in Ref. [1], functional equations for one-loop integrals corresponding to diagrams 
with four external legs can be derived, for example, from the following equation obtained in 
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Refs. [5, 7]: 



Ga {m'(,m2, mi,, mi, m^; {skr}) - {dj/:\^)l''^' {m{,m2,mi,, mi,mi,; {skr}) 



id), 



(djdiA^) li'''\m'^,ml,ml,m^;s23, 834,345,825; 835,824) 
+ {djd2/^5) li'^\ml,ml,ml,ml; 813, 834, 845,815; 835,814) 
+ (Sj^sAs) li'''\ml,ml,ml,ml; 812: 824, 845, 815; 825,814) 
+ {8^84/^5) li'^\ml,ml,ml,ml; 812, 823, 835, 815; 825,813) 
+ (Sj^sAs) li'''\ml,ml,ml,ml; 812, 823, 834, 814; 824,813), 



(1) 



where the operator j"*" shifts the indices i^j^i^j + 1, G4 is the Gram determinant, and A5 is the 
modified Cayley determinant, defined as 
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(3) 



The integral I^\{m'^^, {s^r}) corresponds to a diagram with five external legs and the integrals 
I^\{m^},{skr}) in Eq. (1) are defined as 



(12 2 1 \ 
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(4) 



The diagram corresponding to the integral 
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mi,m2,m^,m 
P23 



q-P3 



\; {^ij}) is presented in Fig. 1. 
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q-P4 



Figure 1: Diagram corresponding to the integral (mj, m^, mg, m^; S12, ■S23, S34, S14; 524,513). 

In what follows, we will use the following short-hand notation for integrals needed in the 
calculation of the one- loop radiative corrections to the process e^e~ — > e~^e~, the so-called 
Bhabha scattering [8], and heavy-quark production: 

B{s,t) = l['^\o,m'^ ,0,m'^; rr? ,rn^ ,rn^ ,vr?; 8,t), 
D2(s,i) = jf^ (0,0,0, m^; 0,0, m'^,m'^; t,s). (5) 
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The diagrams corresponding to these integrals are depicted in Fig. 2. 




Figure 2: Diagrams corresponding to the integrals D2{t, s) and B{s, t). Lines with a zero internal 
mass rui = (for internal lines) or a zero virtuality Sij = (for external lines) are shown wavy. 
Solid lines have a non-zero internal mass or a non-zero virtuality. 



3 Functional equations for the integrals B{s,t) and D2{t, s) 

In this section, we present in detail the derivation of functional equations for the scalar in- 
tegrals B(s,t) and D2{t,s). As was proposed in Ref. [1], one can obtain a functional equa- 
tion for the integral l4^\{'m'j}, {s^r}) by eliminating terms with I^'^\{m'j}, {s^r}) from Eq. (1) 
through an appropriate choice of kinematical variables. The integral -^5*^^ ({m-j}, {sfcr}) depends 
on 15 kinematical variables, while the integral I^\{m'j}, {skr}) depends on 10 variables. There- 
fore, to obtain a functional equation for the integral l['^\{m'j}, {skr}) with all 10 kinematical 
variables arbitrary, we can impose conditions on some 5 variables. To eliminate terms with 
I^'^\{m'j}, {skr}) from Eq. (1), two equations must be fulfilled: 

G4 = 0, djA^ = 0, (6) 

thus fixing two kinematical variables. There are several options to use the remaining three 
variables. One option is to set these variables to some particular values and obtain a func- 
tional equation connecting the integral of interest with integrals which are easy to evaluate 
and/or integrals similar to the original one but with different kinematics. Another option to 
choose variables is to reduce the number of terms in the functional equation by requiring some 
derivatives didjA^ to be zero. Also, one can use some combination of these two options. The 
derivation of the numerous functional equations corresponding to the different options and their 
combinations described above was done on a computer. We describe the most useful functional 
equations below. 

To be definite, let us assume that our integral of interest is the last integral on the right-hand 
side of Eq. (1). If we set in Eq. (1) mi = rus = 0, = ml = m?, S12 = S23 = S34 = S14 = m?, 
and S24 = s, S13 = t, then the last integral on the right-hand side of this equation corresponds 
to our integral B{s,t). 

Setting m| = and choosing different particular values of the remaining two kinematical 
variables and/or requiring some second derivatives didjA^ to be zero, one can get from Eq. (1) 
rather different functional equations for the integral B{s, t). For some specific choice of kinemat- 
ical variables, also the integral D2{s, t) appears in the functional equation. By imposing different 



3 



conditions, wc try to find equations connecting the integrals B{s,t) and D2{s,t) with simpler 
integrals, for example, with integrals having more massless propagators and simpler external 
kinematics. We would like to note that, in the present investigation, we always set m| = in 
order to avoid the appearance of integrals with three propagators having nonzero mass in the 
functional equations. 

Substituting j = 2, ml = = = 0, si2 = S23 = ■534 = S14 = rri^, S24 = s, and S13 = t 
into Eq. (1) and choosing S15, S25, S35, and S45 from the conditions 

G4 = 0, d2A5 = 0, ai52A5 = 0, 52a3A5=0, (7) 

we arrive at the following equation: 

2 2 

B{s,t) = ^{1 + a+)D2{t,m^a+) + ^{1 + a-)D2{t,m^a-), (8) 

where 



Thus, we have a relation connecting the integral -B(s, t) with an integral having only one massive 
propagator, i.e. with the integral D2{s,t). It turns out that the integral D2{s,t) in Eq. (8) 
satisfies the following functional equation: 

D2it,s) = ^D2 (.,!^) +^ If (o, 0,0,0; 0,0,0,0; , (10) 

which can be obtained from Eq. (1) by setting j = m\ = m\ = m| = m| = s\2 = S23 = 0, 
S34 = ■S14 = ml = m^, and S24 = s, S13 = t and imposing the conditions 

G4 = asAs = ai^sAs = dsd^As = 0. (11) 

The last integral in Eq. (10) corresponds to the box integral with all propagators massless and 
the squares of all external momenta equal to zero. By using Eq. (10) and taking into account 
the relation a+a- = 1, one can write the integral in Eq. (8) with argument a+ as 

D2{t, m'^a+) = a- D2{t, m^a-) + (1 - a_) if (0, 0, 0, 0; 0, 0, 0, 0; s- W, t) . (12) 

Substituting this relation into Eq. (8), gives 

S(s,0 = (l-^.)^2(t,m2a_) + /3, if (0,0,0,0; 0,0,0,0; s-4m^,t). (13) 

We illustrate this relation in Fig. 3. Prom Eq. (13), the master integral from heavy-quark 
production is found to be 

D,(s.t) = t±fB i^i^.s) + /f (0.0,0,0; 0,0.0,0; ^^^..) . (14) 

Here, we would like to remark that, for e = (4— d)/2 0, the integral 1^2 has a pole proportional 
to while the leading singularity of the integral B(s,t) is 1/e. The leading 1/e^ singularity 
on the right-hand side comes from the massless integral I4. 

Analytic formulae for both integrals on the right-hand side of Eq. (14) are given in the next 
sections. 
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Figure 3: A schematic depiction of Eq. (13). Wavy lines correspond to massless scalar propaga- 
tors and solid lines to massive propagators. 



4 Analytic result for the integral B{s,t) 

Using the method of dimensional recurrences [5, 6], the following hypergeometric representation 
for the integral B{s,t) was obtained in Ref. [9]: 
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and l2^^ are the one-loop propagator type integrals 
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Here, the Appell hypergeometric functions are 
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The function (j){x,y) is 



(x, y) = fJi^.'o 



1; 1; 



d-3 . d-3 

d-1 . d-2 . _. 

2 ■ 2 ' ' 



00 00 / d-3 ^ / d-3 \ 
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(15) 
(16) 



(17) 



(18) 



(19) 



where F^.''^.^ is the Kampe de Feriet function [10]. The Appell function Fi admits a one- 
fold integral representation (see Appendix). The functions F2 and cp admit two- fold integral 
representations, and this is the reason why their e expansions are problematic. 
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We discovered that both functions can be represented in terms of the Gauss hypergeometric 
function 2F1 and Appell function F3, defined as [10]: 



^ (7)m+n rn\ n\ 



m,n=0 



\x\ < 1, \y\ < 1. 



(20) 



We found two methods to obtain such a representation for the (f) function. The first method 
is as follows. We write the function 4>{x,y) as 

r=0 \ 2 )r \ 2 )r ^ 



1 ^^-3 I », . 
2 ^ ' ' 



(21) 



and then perform an analytic continuation of 2^1 transforming it to two functions 2^1 with 
argument l/y. Thus, we obtain two terms. One of these terms is just the Gauss function 2-F1 
and another one is the Horn function 



Ho 



d-5 d-2 1 

-L 1 J- 1 -L 1 ~ 5 5 



; -^5 



(22) 



defined as [11] 



H2[oi,(3,^,^ ,6,x,y) = 2^ — — ' 0-FI>1' Fl'bl < 1- 



m,n=0 



ml n! 



(23) 

The Horn function (22) can be expressed in terms of 2-^1 functions and the Appell function Fs 
using Eq. (65) on p. 295 in Ref. [12]. Combining all terms, we arrive at the following result: 



Hx,y) = - 
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TT (1 - y)(-x)^^ ^ ^ 



' ' ' 2' ' 2 'x' 1-y 
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(24) 



The Appell function F3 for this specific set of parameters can be written as a one-fold integral: 



r (9^) x{y - 1) (1 - ^ ^'■csin y/^^ 



' ' 2' ' 2 'x'l-yy r(|)r(3-f) 70 {l-x-v)Vl-y-v 



dv. (25) 



This formula can be obtained from the integral representation of the F3 function [11] (see also 
Appendix). The function 0(— m^z, 0) from Eq. (15) reads: 



-m'^z,e) = {d-3)r 



d- 5 



V^T (2 - I) 



(-9)— 2F1 



1, ^ ; 4iv?u 

{t + u){s + u)_ 



2\/w?stu 



t{t + u) 



d-5 
2 



arcsm 



w 
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(26) 
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where 



and 



(pi{s,t) 



I [1 — vj 2 arcsin 



4m2 



4:{u+t) 



dv, 



(27) 



w = 16m^ - Ws - ts. (28) 

From this expression, it follows that any eoeffieient in the e expansion of the (j) function can be 
expressed in terms of one-fold integrals. Several terms of the e expansion of the function (p are 
given in Ref. [13]. 

We present also another method to represent the function (f){x,y) in terms of the Appell 
function F3 and the Gauss functions 2-^1 • For the function (f){x,y), one can write the following 
integral representation [13]: 

1 
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(29) 



From this integral representation, one can derive a differential equation. Differentiating both 
sides of Eq. (29) w.r.t. x, using the following formula for the derivative of the Gauss function. 
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(30) 



z{l-z) 

and after some simplification of the resulting integrand, one obtains the following equation [9]: 
d(f){x,y) 



2{y - x)x- 



dx 



= [y - {d - 3){y - x)](f){x,y) 



(31) 
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This first-order differential equation can be easily solved yielding: 
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(32) 



As it happens in the previous case, instead of the Kampe de Feriet function, we obtain the 
more familiar Appell function F3, which, for the above parameters, admits the one- fold integral 
representation: 
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This relation is obtained from the integral representation given in the Appendix. By using 
Eqs. (32) and (33), we obtain the following expression for the function (f){—m^z, 6) from Eq. (15): 
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where 
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(35) 



9 is defined in Eq. (16) and w in Eq. (28). 

Thus, we obtained two rather different hypergeometric representations for the (j) function. 
In both cases, the hypergeometric functions admit one-fold integral representations, so that all 
the coefficients in the e expansion may be expressed only in terms of one-fold integrals. 

The Appell function F2 from Eq. (15) can also be expressed in terms of the Appell function 
F3 and the Gauss function 2-F1. To obtain such a relation, we use the formula for the analytic 
continuation of the Appell function F3 from Ref. [11] and obtain: 
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The Appell function F3 with this particular set of parameters also can be expressed in terms of 
the one-fold integral 
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Therefore, the Appell function F2 from Eq. (15) reads: 
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(39) 



Thus, we found that the Appell function F2, the function cp and, therefore, also the integral 
B{s,t) are expressible in terms of hypergeometric functions admitting one- fold integral represen- 
tations. Such a representation of B{s,t) would be convenient for obtaining higher-order terms 



8 



in the e expansion of this integrah The first terms in the e expansion of the integral B(s, t) were 
obtained in Rcf. [14]. 

In Ref. [9], a one-fold integral representation for a scalar box integral with arbitrary masses, 
external momenta and space-time dimension d was presented. To obtain the formula for B(s, t) 
directly from such a representation by setting masses and scalar invariants to their specific values 
would be not so easy because the appropriate analytic formulae are rather lengthy, and also the 
analytic continuations needed for the hypergeometric functions are rather nontrivial. Our result 
can be considered as a confirmation that the one-fold integral representations for box integrals 
with physically relevant kinematics do exist. 



5 Imaginary parts and spectral representation for the integral 

The integral with all internal lines massless and external legs on shell can be calculated 
analytically. Using the method of dimensional recurrences [6], we obtain the following relation 
for the last integral in Eq. (13), assuming |s -|- t| < \s\ and |s -|- i| < |i|: 
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(40) 



This formula is in agreement with the result obtained in Ref. [15] (see also Refs. [16, 17]). Thus, 
Eqs. (14), (15), and (40) provide us with a hypergeometric representation for the integral D2- 
The £ expansions of the real and imaginary parts of this integral through order were given 
in Ref. [18]. We expect that, with our hypergeometric representation, one can derive a shorter 
result for the term in the expansion of the integral D2 than that given in Ref. [18]. 

The functional equations in Eqs. (10) and (13) can be used for finding imaginary parts of 
the integrals in some kinematical regions. As one can see from Eq. (10), if s > m? and t < 0, 
the integral l[^\o,0,0,m^; 0,0,m? ,m?; s,t) has an imaginary part that arises only from the 
integral I^'' with all propagators massless, which can be easily found from Eq. (40), so that 
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In a similar fashion, if s > 4m? and t < 0, the imaginary part of the integral B{s,t) can be 
obtained from Eq. (13). In this case, the imaginary part originates only from the second integral 
on the right-hand side of Eq. (13). Again one can use Eq. (40) to find: 
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One can write D2{s,t) in the fixed-t spectral representation 
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Substituting Eq. (41) into Eq. (43) leads to the following integral representation for D2{t,s): 
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To use D2{t,s) in calculations of heavy-quark production, one needs to know it at t > and 
s < 0. To perform the analytic continuation of D2{t,s) into the region t > 0, one can use 
Eq. (44). Using a formula for the analytic continuation of the Gauss hypergeometric function 
[11] and introducing the new integration variable 
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leads to the following expression: 
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where 



2\2 



(7 = 1 + 



(s — m ) 
'st 



(46) 



(47) 



and d is defined in Eq. (16). Substituting the expression for in this formula, we obtain for 
the imaginary part: 



Im L>2(i,s) 



-7r2(d-3)t2- 
s2d-iT (^) I d-4 



2^1 



d-2 .(^ 
2 ' J 



s + m 
d-3 



F, 1,1,^ 



1 d- 1 



' '2' 2 



;(7,e 



s < 0, t > 0. 



(48) 



Applying the Eulcr transformation to the hypergeometric functions 2F1 and Fi (see, for example, 
Ref. [11]), we find agreement with the expression for the imaginary part obtained in this region 
from the analytic result for the integral B(s,t) given in Eq. (15). 
For the real part, one can write the following representation: 



I ^ • ■ 

d-2 . ^ 
2 1 



+ coty ImD2(t,s) + $(i,s), 

s < 0, t > 0, 
(49) 
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where 



, 2^-'^ (d - ?,)is + m^) } dv v^-^il-v)^ ^ 
^\t,s) = ^ — ^ / / — 2F1 



1 — ^ • 

' 2',, 

^ 2 ' . 



(50) 



Equations (48) and (49) are convenient for the e expansion. Appropriate formulae for the e 
expansion of the Gauss hypergeometric function can be taken from Refs. [19, 20]. To obtain 
the £ expansion of the Appell hypergeometric function Fi, one can use a one-fold integral 
representation (see Appendix). The first two terms in Eq. (49) give singular contributions 
proportional to 1/e^ and 1/e. The last term in Eq. (49) is regular in e. Its £ expansion can be 
derived from integral representation of Eq. (50). 

In the case when t > and s < 0, the imaginary part of B(s, t) can be found from Eq. (15) 
and reads: 



7ri 2^-"' _ /d - 3 1 d - 1 



Im B(s,t) = t^t-ttFi , 1, -; ; 1 ^ , 1 ^ . (51) 

The imaginary part of the integral D2{s,t) for s > and t < reads: 



Imi)2(.,i)=sin^4'^)(0,0; - s) j—^^-^-— 2F, 



1,1-2;-, , St 



2 



1 + 



(t + m^) fd-3 1 d-l st s 



m^{t-m^)^^^\ 2 '-^'2' 2 '^+(^-^2)2'^ 4m2 J j ' ^^^^ 
The first few terms of the e expansion of the Appell function Fi can be found in Refs. [9, 21]. 

6 Functional equations for D2 with one quark leg off shell 

For the computation of the radiative corrections to heavy-quark production in the NNLO ap- 
proximation, one needs to know the integral D2 with one quark leg off shell. This is a rather 

nontrivial task, which will be considered in a forthcoming publication. In this section, we just 
want to outline the strategy of such calculations. Similar to the case of D2 with all legs on shell, 
one can derive the functional equations: 

/f^(0,0,0,m2;0,0,S34,m2; s,t) 

= —lf 0,0,0,m^0,^ ^ L^su,w?.—,t 

s \ s s 

s — m? jd) ( „ „ „ ivr? — S34)(m2 — s) „ „ „ (vr? — s? , , „, 

+ if 0, 0, 0, 0; ^ ^, 0, 0, 0; ^ ^-,t , 53 

s \ s s ' 

/f)(0,0,0,m2; 0,0,S34,m2; s,t) 



S34 \ S34 S34 S34 



+ fo,0,0,0; ^,^^^^^^^,0,0; i^' ' ss,){m' - s) \ _ 

S34 \ S34 S34 S34 J 
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Equations (53) and (54) were obtained from Eq. (1) by setting 

2 2 2 2 2 2 

rUi = 1712 = "^3 = "T-S = ■^12 = S23 = 0, 777,4 = ■^14 = '^^ , 

G4 = 55A5 = Si^gAs = 0, (55) 
in both cases and additionally 

aa^sAs = 0, (56) 

in Eq. (53) and 

^s^sAs = 0, (57) 

in Eq. (54). The above equations are rather similar to the functional equations for the integral 
Z?2 with all legs on shell. In some kinematical regions, the imaginary part of the integral 
D2{ss4,s,t) arises from integrals with all propagators massless. In Eq. (53), the integral with 
massless propagators can be expressed in terms of three Gauss hypergeometric functions [15]. 
The analytic expression for the massless integral in Eq. (54) for arbitrary d is not known at 
present. The first terms in the e expansion of both integrals may be found in Refs. [15, 16]. 
We expect that, similar to the D2 integral with all legs on shell, the integral D2 with one leg 
off shell can be represented in terms of hypergeometric functions admitting a one-fold integral 
representation. 



7 Conclusions 

The usefulness of functional equations turns out to be threefold. First, we obtain a hypergeomet- 
ric representation of the integral needed for NLO calculations of heavy-quark production from 
the result for the integral from Bhabha scattering. Second, since the term in the expansion 
of the integral D2 is known [18], one can use it to obtain the term for Bhabha scattering. 
Third, functional equations provide a simple method to obtain imaginary parts of integrals. For 
some kinematic regions, the imaginary parts of integrals with nonzero internal masses can be 
related to the integrals with all lines massless. 

It is also important that the functional equations provide a tool for performing analytic 
continuations of the considered integrals. As was already observed in Ref. [1] and now also in 
this paper, functional equations suitable for analytic continuation express the considered integral 
in terms of the same integral with transformed arguments that has no imaginary part plus simpler 
integrals (usually with massless lines) giving the imaginary part of the integral. Such analytic 
continuation is achieved by solving algebraic equations, so that the explicit analytic form of the 
integral is not needed. 
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9 Appendix 

In this appendix, we present formulae which were used in the derivation of some equations in 
the main text. 
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The integral representation of the Appell hypergeometric function Fi reads: 

1 

-Fi (a, /?,/?', 7; x,y) = / -J— -grdu, 



Re a > 0, Re (7 - a) > 0. (58) 

The integral representation of the Appell function F3 used in the derivation of the one-fold 
integral representations of the (f)i functions reads: 

F2{a,a',P,P',r,x,y) = —^^^^-— ^^/^ '''^~''\[ \ f, 

r(7 - /3)r(/3) 7o {l-x + ux}'^ 

This integral representation follows from Eq. (20) in Ref. [22] . 
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